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Abstract 



The order a s corrections to the Wilson coefficients of the dipole operators (O7, 0$) at the 
matching scale /z = my/ are a crucial ingredient for a complete next-to-leading logarithmic 
calculation of the branching ratio for b —>■ 37. Given the phenomenological relevance and the 
fact that this two- loop calculation has been done so far only by one group [|TJ, we present a 
detailed re-calculation using a different method. Our results are in complete agreement with 
those in ref. |l|. 
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1 Introduction 



By definition, rare B meson decays only arise at the one loop level in the standard model (SM). 
Therefore these decays are particulary sensitive to effects from new physics. Among these 
decays, the inclusive modes like B — > X s ^ are particulary interesting, because no specific model 
is needed to describe the final hadronic state in contrast to the exclusive decay modes. Indeed, 
heavy quark effective theory tells us that the decay width T(B — > X s 7) is well approximated by 
the partonic decay rate T(b — > X s / y) which can be analyzed in renormalization group improved 
perturbation theory. The class of non-perturbative effects which scales like is expected 

to be well below 10% Q. This numerical statement is supposed to hold also for the recently 
discovered non-perturbative contributions which scale like l/m 2 c ||. 

Up to recently, only the leading logarithmic (LL) perturbative QCD corrections were cal- 
culated systematically M. The error of these calculations is dominated by a large renormal- 
ization scale dependence at the ±25% level. The measured branching ratio BR(B — > X^) = 
(2.32 ± 0.67) x 10~ 4 reported in 1995 by the CLEO collaboration |J overlaps with the esti- 
mates based on leading logarithmic calculations (or with some next-to-leading effects partially 
included) and the experimental and theoretical errors are comparable |, 0, B B S3 Q- How " 



ever, in view of the expected increase in the experimental precision in the near future, it became 
clear that a systematic inclusion of the next-to-leading logarithmic (NLL) corrections becomes 
necessary [[J. This ambitious NLL enterprise was recently completed; combining the results of 
different groups |], ||, [10], [12], [13], |14| , the first complete theoretical prediction to NLL pecision 



for the b -> X s + ^ branching ratio was presented in JTJ: BR(B -> X„i) = (3.28 ±0.33) x 10" 4 . 
This prediction is still in agreement with the CLEO measurement at the 2cr-level. The the- 
oretical error is twice smaller than in the leading logarithmic prediction. So the inclusive 
B — > X s ± 7 mode will provide an interesting test of the SM and its extensions when also more 
precise experimental data will be available. 

Before discussing in some more detail the principle steps leading to a next-to- leading result 
for b — > Xg'-f, we briefly have to recall the formalism. We use the framework of an effective low- 
energy theory with five quarks, obtained by integrating out the top quark and the W^-boson. 
The effective Hamiltonian relevant for b —>■ s'j and b — > sg reads 

AC 8 

H eff {b^s 1 ) = —£\ t Y,C i {Li)O i {v) , (1.1) 
V2 i=1 

where Oj(/z) are the relevant operators, Cj(/i) are the corresponding Wilson coefficients, which 
contain the complete top- and W- mass dependence, and \ t = V t bV t * s with being the CKM 
matrix elements 2 . Neglecting operators with dimension > 6 which are suppressed by higher 
powers of l/mvK/t-factors and using the equations of motion for the operators, one arrives at 
the following basis f\ of dimension 6 operators ]15 



0\ = {c L ^b La ) (s Lq 7 m c L/3 ) , 
2 = {c La l^h La ) (s L/3 7 M c L/3 ) , 
O-z = (s La ^b La ) {ulpIhUlp) ± ••• + (bLpJvhp 



2 The CKM dependence globally factorizes, because we work in the approximation X u = 0. 
3 In |14| another basis was used. We comment on this in the summary. 
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7 = (e/lQTT 2 )s a a^(m b (fi)R + m s (fi)L)b a , 

Os = {g s /\QT: 2 )s a a^{mb{ii)R + m 8 {ii)L){\%/2)bpG% . 

(1.2) 

In the dipole type operators O7 and Os, e and F^ u (g s and G^,) denote the electromagnetic 
(strong) coupling constant and field strength tensor, respectively. 

It is well-known that the QCD corrections enhance the b — > 37 decay rate by more than a fac- 
tor of two; these QCD effects can be attributed to logarithms of the form a;™(ra&) log m (m^/M), 
where M = m t or M = my/ and m < n (with n = 0,1,2,...). Working to NLL preci- 
sion means, that one is resumming all the terms of the form a"(m(,) ln n (m&/M), as well as 
<^s{ m b) i a ™{ m b) hi n (m^/M)). This is achieved by performing the following 3 steps: 

Step 1 One has to match the full standard model theory with the effective theory at the 
scale \i = fiwu where fiwt denotes a scale of order mw or m t . At this scale, the matrix 
elements of the operators in the effective theory lead to the same logarithms as the full 
theory calculation. Consequently, the Wilson coefficients Ci(jiwt) only pick up small 
QCD corrections, which can be calculated in fixed order perturbation theory. In the NLL 
program, the matching has to be worked out at the 0(a s ) level. 

Step 2 Then one performs the evolution of these Wilson coefficients from fi = fiwt down to 
fi = fib, where fib is of the order of m,b- As the matrix elements of the operators evaluated 
at the low scale //& are free of large logarithms, the latter are contained in resummed form 
in the Wilson coefficients. For a NLL calculation, this RGE step has to be performed 
using the anomalous dimension matrix up to order a 2 s . 

Step 3 The corrections to the matrix elements of the operators (s7|Oj(/i)|6) at the scale = 
have to be calculated to order a s precision. 

The most difficult part in Step 1 is the two- loop (or order a s ) matching of the dipole operators, 
which has been worked out by Adel and Yao U some time ago. Step 3 basically consists of 
Bremsstrahlung corrections and virtual corrections. The Bremsstrahlung corrections, together 
with some virtual corrections needed to cancel infrared singularities, have been worked out 
by Ali and Greub || ID|; later, this part was confirmed and extended by Recently, a 



complete analysis of the virtual corrections (up to the contributions of the 4 Fermi operators 
with very small coefficients) were presented by Greub, Hurth and Wyler [13]. The main result 



of the latter analysis consists in a drastic reduction of the renormalization scale uncertainty 
from about ±25% to about ±6%. Moreover, the central value was shifted outside the la bound 
of the CLEO measurement. However, at that time, the essential coefficient Cr(fib) was only 
known to leading-log precision. It was therefore unclear, how much the overall normalization 
will be changed, when using the NLL value for C 7 (fj,b). Very recently, the order a 2 anomalous 
matrix (step 2) has been completely worked out by Chetyrkin, Misiak and Miinz [|14j]. Using 
the matching result of Adel and Yao, these authors got the next-to- leading result for C7 (//&). 
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Numerically, the LL and the NLL value for CV(//&) are rather similar; the NLL corrections to 
the Wilson coefficient C 7 (//&) lead to a change of the b — > X s ^ decay rate which does not exceed 
6% : The new contributions can be split into a part which is due to the order a s corrections 



to the matching (Step 1) and into a part stemming from the improved anomalous dimension 
matrix (Step 2). While individually these two parts are not so small (in the NDR scheme, which 
was used in |14|]), they almost cancel when combined as illustrated in |[l4j]. This shows that all 



the different pieces are numerically equally important. However, strictly speaking the relative 
importance of different NLO-corrections at the scale /i = fit>, namely the order a s corrections 
to the matrix elements of the operators (Step 3) and the improved Wilsoncoefficients Cj (Step 
1+2), is a renormalization-scheme dependent issue; so we stress that the discussion above was 
done within the naive dimensional regularization scheme (NDR). 

Each of the three steps implies rather involved computations: The calculation of the matrix 
elements (Step 3) involves two-loop diagrams where the full charm mass dependence has to be 
taken into account. Also the matching calculation (Step 1) involves two-loop diagrams both in 
the full and in the effective theory. Finally, the extraction of some of the elements in the 0{a 2 s ) 
anomalous dimension matrix involves three-loop diagrams. Given the fact, that it took a rather 
long time until the leading logarithmic calculations performed by different groups converged 
to a common answer, it is certainly desirable that all three steps mentioned above should be 
repeated by other independent groups, and, may-be using other methods. 

Making a step into this direction, we present in this paper a re-calculation of the two-loop 
matching of the dipole operators O7 and 0%. We extracted the 0(a s ) contributions of the 
corresponding Wilson coefficents C7 and C% by calculating the on-shell processes b — > sj and 
b — > sg in both versions of the theory up to order a s . We worked out the two-loop integrals by 



using the Heavy Mass Expansion method fL6fl , which we describe in some detail in section 2.4. 

The rest of the paper is organized as follows. In section 2 we make some preparations 
for the two-loop calculations. We first explain how to extract the order a s corrections to the 
Wilson coefficients C-j{fiwt) and C%(nwt) in principle. Then, in various subsections we discuss 
and illustrate the technical methods used. Sections 3, 4 and 5 are devoted to the computation 
of Cni^wt)'- in section 3 we calculate QCD corrections to b — > S7 in the full theory together 
with the corresponding counterterm contributions, while in section 4 the same is done in the 
effective theory. Comparing the results from section 3 and section 4, we extract Cn(nwt) in 
section 5. Similarly, sections 6, 7 and 8 are devoted to the computation of Csi(fiwt)'- In section 
6 we calculate QCD corrections to b — > sg in the full theory together with the corresponding 
counterterm contributions, while in section 7 the same is done in the effective theory. Comparing 
the results from section 6 and section 7, we extract Csi(fJLwt) in section 8. Finally, we give a 
brief summary in section 9. 

2 Preparations for the two-loop calculations 



2.1 Strategy for extracting C71 and Q 
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Let M denote the (on-shell) b — > 37 matrix element calculated in the effective theory. M can 
be written in the form 

M = E^)(°W) » (Ofo)) = (*y\OiOA)\b) . (2.1) 



3 



To keep the notation simpler, we denote the matching scale by fi instead of jj, wt . Making use 
of the a s expansion for Cj(//) and Oj(/x) 

C^) = CM + ^C a ^) + ... , (O i (ji)) = (O i (p)) + -£(O i ( f i))i + ... , (2.2) 

we get the corresponding expansion for M in the form 

M = CM {Oi(ji)) Q + / {CM (Oi<J*))i + CM (Oi(ji)) Q ) + ■■■ ■ (2.3) 

On the other hand, let M denote the b — > S7 matrix element evaluated in the full theory 
after discarding power supressed terms of order l/mL t ; M has the expansion 

M = M + ^ Mi + . . . . (2.4) 

47T 

Requiring M = M and taking the coefficient of a\, we get the 0(a s ) matching condition 

M 1 = CM(Oi(»))i J rCM(Oi{v))v • (2-5) 

All coefficients in eq. ( |2.5|) are known |L7| pL8|| , except [] C71 and C 8 i- As C 8 i comes together 
with (sy\0$(fjL)\b)o, which is zero, eq. ( |2.5|) has only one unknown, viz. C71, i.e., just what we 
want to extract. 

The discussion for the extraction of Cgi goes exactly along the same lines, using the process 
b — > sg instead of b — > 37. 

A general remark is in order here. One could also match off-shell Greens functions instead 
of on-shell matrix elements. However, in this case one is not allowed to work in the operator 
basis given in eq. (|1 . 2|) , because one has used the equations of motion for the operators to get 
this 8 dimensional basis. This Hamiltonian therefore only reproduces on-shell matrix elements 
correctly ||19|| . As we would have to work in the off-shell basis when matching Greens functions, 
we preferred to do on-shell matching. There is of course a price to pay: The on-shell processes 
b — > 57 and b — > sg are plagued with infrared singularities, which have to be treated carefully. 
However, as we will see later, this is not a real problem. 

2.2 Technical details 

We work in d = 4 — 2e dimensions; in the full theory we use anticommuting 75, which should 
not be a problem, because there are no closed fermion loops involved. We also use this naive 
dimensional regularization scheme (NDR) in the effective theory. The calculations are done 
in the 't Hooft-Feynman gauge (electroweak sector) and the gluon propagator is taken in the 
Feynman gauge. To avoid Euler 7^ terms and ln(47r) factors in our expressions, we introduce 
the renormalization scale in the form fi 2 exp(7 S )/(47r) (MS subtraction then corresponds to 
subtracting the poles in e). Besides the ultraviolet singularities also the infrared singularitites 
are dimensionally regularized. As we could clearly separate infrared and ultraviolet singularities, 
we labeled the infrared poles by the index ir (e.g., l/e^). We put m s = 0, except in situations 
where mass singularities appear, i.e., we treat m s as a regulator of these singularities. We 
work in the approximation A n = 0. To keep the formulae more compact, we put immediately 
Q u = 2/3 (Qd = —1/3) for up-type (down-type) quark charges. For the same reason we also 
immediately insert the numerical values for the color factors in the b —>■ sg case. 

4 Of course C71 and Csi are also known from Adel and Yao M, but this is what we want to check. 
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2.3 Reducing the number of diagrams 

For reasons of gauge invariance, we know that the final result for the b — ► 37 matrix element 
can be written in the form 

M(b — > 37) = F(masses, couplings) (s^\0-j\b) tr ee ■ (2.6) 

For m s = 0, the quantity (s^\0 7 \b) tr ee is given by 

16 7r 2 r i 

(s7|0 7 |6) iree = 2m b u(p') $ 4 Ru(p) = u(p') 2m\i L - 4m 6 (pe) i? u(p) , (2.7) 

e L J 

where tt(p') and tt(p) are the Dirac spinors for the s and the 6 quarks, respectively, and q (e) 
the momentum (polarization vector) of the photon. In the last step we used q = (p — p') and 
qe = , where p (p') is the momentum of the b- (s-) quark. When calculating a given Feynman 
diagram, it is sufficient to work out only the term proportional to (pe) R. After adding all the 
diagrams, the full answer can be reconstructed by means of eq. ( |2.7| ) . This reduces the number 
of diagrams; e.g., when calculating the 0(a s ) corrections for b — > 57 in the full theory, "only" 
the graphs in Fig. 2 have a non-zero projection on the term (pe)R. 

A similar projection for the process b — > sg (with obvious changes) can also be obtained. 

2.4 Method for calculating of the two-loop diagrams 

To extract C71 and Cgi various one- and two- loop diagrams have to be calculated in both 
versions (full/effective) of the theory. As the one-loop diagrams are straightforwardly obtained 
by conventional techniques, we directly move to the two-loop diagrams. When working out 
b — > S7 and b — > sg in the effective theory at the matching scale yLwu the only two- loop 
contributions leading to terms of order a s are those associated with the operator Oi- For the 



b — > 57 case, these terms have been obtained in |fL3| . We anticipate, that in the corresponding 
full theory calculation a term appears which can be identified with the O2 contribution in the 
effective theory. Consequently, the O2 contribution is not needed explicitly for extracting C71 
and C%\. 

Therefore, we directly discuss the calculation of the two-loop contributions in the full theory. 
In order to match dimension 6 operators, it is sufficient to extract the terms of order m;, ^ 
(M = mw, Trit) from the full-theory matrix elements for b — > 57 and b — > sg (term supressed by 
additional powers of rrib/M correspond to higher dimensional operators in the effective theory). 
A systematic expansion of the matrix elements in inverse powers of M can naturally be obtained 
by using the well-known Heavy Mass Expansion (HME). In our context we use this HME only 
as a method for working out the dimensionally regularized two-loop Feynman graphs (and not 
to get directly renormalized quantities). The theory of asymptotic expansions of Feynman 
diagrams is already a textbook matter [BD| 0. Therefore, we only recall those properties of the 
HME, which are of practical importance for our calculation (for the mathematical foundations 

5 The idea of deriving operator product expansions using subtractions of leading asymptotics goes back to 
Zirnmermami [El|. Later this idea was systematically developed within the BPHZ scheme p2f. The simple 



explicit formulae for asymptotic expansion within dimensional regularisation like (2.8) have been systematically 
derived in IfLql. 
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of this method we refer to the literature [|T6] ): 



Suppose that all the masses of a given Feynman diagram T can be divided into a set of large 
M = {Mi, M2, . . .} and small m = {mi, rri2, ■ ■ ■} masses and assume that all external momenta 
Q = {Qi-, Q21 ■ ■ •} are small compared to the scale of the large masses M; then the statement is 
that the dimensionally regularised (unrenormalized) Feynman integral associated with the 
Feynman diagram V can we written as 

Fr M ~°°E f r/ 7 o^^F 7 (f,mM) , (2.8) 

7 

where the sum is performed over all subgraphs 7 of V which fulfill the following two conditions 
simultaneously: 

• 7 contains all lines with heavy masses (M) and 

• 7 consists of connectivity components that are one-particle-irreducible with respect to lines 
with small masses (m). 



Here some clarifying remarks are in order: 

• The operator T performs a Taylor expansion in the variables q 1 and m 7 , where m 7 denotes 
the set of light masses in 7 and q 1 denotes the set of all external momenta with respect to the 
subgraph 7; to be more specific, an external momentum with respect to the subgraph 7 can 
be an internal momentum with respect to the full graph V. -Fr/7 denotes the Feynman integral 
corresponding to the reduced graph r/7. Note that the operator T is understood to act directly 
on the integrand of the subgraph 7. The decomposition of the original, say /-loop-diagram V 
into the subdiagram 7 and the diagram r/7 is achieved in the corresponding Feynman integral 
by factorizing the product of scalar propagators as Ilr = rir/ 7 n 7 such that 

F r / 7 o Tgy^Fj = J dkx ■ ■ -dki n r/7 Tqi^ylLy. (2.9) 

• The full graph T is always a subgraph contributing in the sum 

• It is instructive to look at the special case where all masses are large compared to the external 
momenta in a given diagram T. In this case only the full graph T contributes to the sum J2y 



in (2J3). The complete HME expansion reduces to a naive Taylor expansion in the external 



momenta of the integrand of the Feynman integral: 



M- 



F r T g r F r (q L , M) (2.10) 

• The Taylor operator T introduces additional spurious IR- or UV-divergences in the various 
terms of the sum 2 7 > as we will see i n an explicit example below. It is a nontrivial property of 
the HME that all these artificial divergences cancel after making a summation over all possible 
subgraphs 7. For our calculations this property provides a nontrivial check for the individual 
diagrams, as this cancellation has to happen diagram by diagram. 



Now we illustrate this rather formal description for the diagram in Fig lb, for an internal top 
quark and denote it D top - It belongs to the Setl in Fig 2. The W- and ^-exchange contributions 
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are understood to be added. The corresponding Feynman integral has the following form (the 
Dirac spinors u(p') and u(p) are amputated): 



Dtop 



Diracu 1 Dirac2t 1 



(j) — q + r) 2 — m 2 r 2 [((I + r) 2 — m 2 ) (I 2 — m 2 ) ((I + q) 2 — m 2 ) ] (I + p — g) 2 — m 2 ^ 

In ( p. 11|) the functions Diracu and Diracit are the respective Dirac structures, whose explicit 
form is not important for explaining the principle steps of the expansion. The constant X 
collects all the remaining constant factors like coupling constants and CKM-factors. 

We find two subdiagrams 7 of D top which fulfill the two conditions given below eq. ( |2.8|) : 
The first contribution of the HME corresponding to the subdiagram 71 shown in Fig. lc is 
given by 

D top = X exp(27*e)A47r)-*/ ^ / ^ * (2-12) 
Diracu 1 ^ ( Dirac2t 1 \ 



Dl p = X exp^KW"- — - — - x (2.13) 



(p — q + r) 2 — m 2 r 2 r,p,q \ [((I + r) 2 — m 2 ) (I 2 — m 2 ) ((/ + q) 2 — m 2 ) ] (/ + p — q) 2 — m 2 ^ 

The second contribution is the naive one, 72 = D top (see Fig. Id): 

d d r r d d l 
(27r) d J (2tt 

/ Diracu 1 Dirac2t 1 \ 

P ' 9 \ (P ~ 9 + r ) 2 — m s r2 [((^ + r ) 2 — m t ) (^ 2 — m t) ((^ + ?) 2 — m t) ] (I +p — q) 2 — mly J 
So we end up with 

Aop-^a^+a; • (2.i4) 

The integrals are considerably simplified after the Taylor operation T and can be solved 
analytically after introducing Feynman parametrization. We mention that the Dirac algebra 



has been done with the algebraic program REDUCE |23| and the integrals have been done with 
the symbolic program MAPLE p4 . 



As mentioned above we can discard terms of order 1/M 3 . Simple dimensional arguments 
tell us that we have to perform the Taylor operation T up to second order in the external 
momenta r, p, q in T)\ and also up to second order in p, q in Df op . Restoring all the factors 
which we symbolized by X, and projecting on the term (pe) R, we get 

Dtop = ^0 g CV ^3 (-4m*) (pe)R [d] op + d 2 top ] . (2.15) 

The quantities d\ op and d 2 top are given by (z = {m t /m w ) 2 )\ 



g 9^2 _ g z _|_ 2 

— (1 - 2 \n(m b /m w ) e - 2 ln(m s /m 6 ) e + 4 \n(fi/m w ) e) — — — 

e 6b{z — l) 4 

llz 4 - Uz 3 + 234z 2 - 180z + 24 , 229z 3 + 15z 2 + 7Uz - 53£ 



lnz — — — (2.16) 



108U-1) 4 " 648(z-l) 3 



„ S , A , , . N N -108z In 2 + 6O2 4 - 2582 3 + 4682 2 - 2942 + 24 / u \ 4e 

<p = +— 1 + 4 In (fi m w )e) + — - — M— 

^ e ir 108(2;— l) 4 e \m w J 

9z 2 + IO2 + 2 , 2 1422 4 - 538z 3 + 7532 2 - 2182 + 26 , 

H ; n — hi 2 : — A In 2 

12(z -l) 4 108(2 -l) 4 

z 4 - 402 3 + 27z 2 - 10 z - 2 1 672 3 + 23432 2 - 27662 + 230 

Li(l ) — — — , (2.17) 



18(z-l) 4 v z' 648(2 -l) 3 

where the function S is 

(-542 2 + 48z - 12) In z + II2 4 - Uz 3 + 27 z 2 -38,2 + 14 

S — ; — ; . (2.18) 

108(2 -l) 4 v ; 

The 1/e poles in dj op correspond to spurious ultraviolet singularities produced in the re- 
integration after expanding the subdiagram 7^ The l/e ir poles in d 2 op on the other hand arise 
due to the worsened infrared behaviour induced when expanding the s-quark propagator. As 
we explicitly see, these artifical singularites cancel when adding dj op and d 2 op . 

We now discuss the corresponding diagram D c h arm where the internal top quark is replaced 
by the (light) charm quark. The quantities D\ harra and D^ harm corresponding to the subdiagrams 
71 and 7 2 (see Fig. lc,d) are given by the analogous formulae ( |2.12| ) and ( |2.13| ), where m t is 



replaced by m c and the Taylor operator T r<Ptq in (|2.12| ) is replaced by % }P}qtmc and T p q in ( |2.13|) 
by T Ptq ^ mc . As we are discarding terms of order 1/M 3 , it turns out that only the zeroth order 
term in the m c expansion has to be retained; this amouts to putting m c = in D\ harm and 

charm' 

Moreover, in the charm-case there is a third contribution to the HME which corresponds to 
the subdiagram 73 in Fig. le. The latter consists of the l^/$-line only. As we neglect terms 
of order 1/M 3 , the Taylor expansion of the corresponding Feynman integral just amounts to 
replace the W and $ propagator by ijvn^ and — z/m 2 ^, respectively. As the Feynman integral 
of the $ diagram has an additional factor of order (m c rrib) / from the Yukawa couplings, 
only the four Fermi version of the W exchange diagram effectively contributes to D 3 harm . Stated 
differently, D 3 harm is directly related to the 2 contribution in the effective theory. Of course, 
this is not suprising when keeping in mind how the effective Hamiltonian is constructed. 

To summarize, D c h arm is given by 

Dcharm ~ D chaTm + D charm + D charm ■ (2.19) 

2.5 The matching to leading-log precision 

To establish some lowest order matching functions which are frequently used in the following 
sections and in order to explain an important subtelty in the NLL matching calculation, we 
recall the results of the LL matching: In the full theory the lowest order matrix elements Mo 
for b —>■ sj and b — > sg are obtained by expanding the diagrams shown in Fig. la up to second 
order in the external momenta. The results read in d = 4 — 2e dimensions 

M (b -» 57) = AlG J 1 Xt K m ( S1 \0 7 \b) , M {b^sg) = A%G J. Xt K 80 (sg\0 8 \b) , (2.20) 
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where the functions K 70 and K 80 have an expansion in e of the form 

K 70 = K 700 + e K 701 + e 2 K 702 + . . . , K 8Q = K 80Q + e K 801 + e 2 K 8m + . . . . (2.21) 

On the other hand, the lowest order result M in the effective theory reads (also in d = 4 — 2e 
dimensions) 

M (b - s 7 ) = C 70 (s 7 |0 7 |6) , M (6-^sp) = C 80 HO 8 |6> . (2.22) 

As the matching is understood to be done in 4 dimensions, we get the connections 

C70 = K 700 , Cso = i^soo • (2.23) 

Therefore, in d dimensions Mq and Mo differ by terms of order e. This detail becomes an 
important subtelty when going to higher loop orders; we best explain this by means of an 
example: one type of order a s corrections is given by multiplying the lowest order result by 
ultraviolet singular v^2 factors which are the same for both versions of the theory. In the full 
theory, this leads to finite terms proportional to K 7m (and K 8m )\ the corresponding terms in 
the effective theory are not generated. 

When working out the two-loop integrals corresponding to the diagrams in Figs. 2 and 3 
in the full theory for b — > s 7 (or b — > sg), there are contributions in which the dimensionally 
regularized lowest order result, taken up to first or second order in e, factorizes. As we will 
see later, the infrared singularity stucture is precisely of this form. As we will use the explicit 
expressions for the Inami-Lim [^J functions K 700 , K 800 , K 7m and K 8m at several places, we 
list them here. Using z = (mt/mw) 2 , they read 

z[6z(3z-2) lnz-Qs-l) (8z 2 + 5z - 7)] 
K 700 = C 70 = ^— — (2.24) 

z[6z \nz + (z-l)(z 2 -5z-2)] 

-"■800 — tv80 — g/ _ ^ 4 {Z.ZO) 



K> 



18z{3z - 2) In 2 z + (44z 3 - 314^ 2 + 324^ - 96) In z + 56z 3 - 35z 2 - 56z + 35 



701 



UA(z - l) 4 

+2K 70Q ln(fx/m w ) (2.26) 



-18z In 2 z + (10z 3 - 28z 2 + 108z - 48) In z + 25z 3 - 118^ 2 + 119 z - 26 



S01 



48(z- l) 4 

+2 K 800 In (avW) (2.27) 

3 6 — >• S7 in the full theory 

In section 3.1 we present the results for the dimensionally regularized matrix element M for 
b — > s 7 in the full theory. In section 3.2 we discuss the various counterterm contributions. 
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3.1 Two- loop Feynman diagrams 



As in eq. (|2.4j ), we write the b — > 37 matrix element M in the form M = Mq + ^ Mi. When 
using the "reduction technique" described in section 2.3, the complete list of two-loop diagrams 
contributing to Mi is given in Fig. 2, where the cross stands for the possible locations where 
the photon can be emitted. Note that diagram 5b in Fig. 2 does not contribute in the limit 
m s = 0. We write the result for Mi in the form 



An 



Mi — V R] +2 - R\ +2 - Rl , (3.1) 



where V is an abbreviation for the often occurring quantity 

V=^^^-C F ( S1 \0 7 \b} tree ; C F = ^ . (3.2) 

In eq. ([□]) R\ +2 {R\ +2 ) denotes the sum of the first and second contribution in the Heavy Mass 
Expansion (HME) (see section 2.4) of the dimensionally regularized (unrenormalized) Feynman 
integrals for internal top (charm) quark; R% is the third contribution in the HME, which has to 
be considered only for the light internal quarks, which in our present case is the charm quark 
(X u = 0). According to the HME, R? c is obtained by working out the charm loops using the 
four- Fermi approximation of the H^-propagator. Stated differently, R? c is directly related to the 
order a s contribution of matrix element of the operator 2 , provided the latter is evaluated in 
the NDR scheme; more precisely, 

R 3 C = -R 2 (3.3) 

where R 2 is the quantity defined through the equation 

cy 

(s 1 \0 2 \b) = -^C F (s 1 \0 7 \b) tree R 2 . (3.4) 

47T 

As the same contribution is also present in the effective theory, we will not need to know R% 
explicitly [] in order to extract the order a s corrections in the Wilson coefficient C 7 i(/j,wt)- 
Making use of the various K— functions given in (|2.24|) - (|2.27 ) and denoting r = (m s /m{,) 2 , we 



now give the dimensionally regularized expression for R 1+2 = R] +2 — R\ + 



2 



R 1+2 = -{K 7m + eK 7m )^^\nr + gi ^-+ l -K im ^ 2 r 

e ir e 2 

+2K 700 \nrln(m b /m w ) -2K~ 700 \nr + g 2 ln(m b /mw) + g3 ■ (3.5) 



The first term in eq. (|3.5| ) is due to infrared singularities in the on-shell b — >• 37 amplitude as 
suggested by the notation [] l/e ir . This term is entirely due to those diagrams in set 3 of Fig. 
2 where the photon is radiated from the internal quark or the W (or <3>) boson. The quantities 
gi, g 2 and #3 in eq. ( |3.5| ) can be written as (z = (m t /mw) 2 , Li (a;) = — Jq y ln(l — t) ) 



9\ 



t 

-324z 4 - 45(k 3 + 270z 2 + 72z) \nz + l 12z 5 + 2Uz 4 + 55z 3 - 931z 2 + 5932 - 73 35 

72(z- l) 5 216 

(3.6) 



6 The reader who whishes to see the explicit form for R2 is referred to eq. (2.35) in ref. Jl3[ . 
7 We could separate ultraviolet and infrared poles in our calculation. In the follwing, l/e, r stands for 

infrared (ultraviolet) poles. 
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92 



-21Qz 3 + 162/ - 72z) In z + 44/ + 154/ - 393/ + 274z - 79 7 

36(^- l) 4 4 



(3.7) 



93 



Ll 1 - -) + - l 7T Kgoo + — 7T 2 + 



1 

2' 



2 
3 



2 

27 



1296 



z(8z 3 + 61/ - AOz + 4) 
6(z-l) 4 

- (-4860/ - 18954/ + 11502/ + 648z) In 2 z 

+ (3240/ + 16956/ + 37638/ - 56586/ + 20688z - 2496 - 216 /(/ - z 2 )) In z 

+ (-1442/ - 55910/ + 109651/ - 69271/ + 20999,2 - 4027) 

+ (60/ - 228/ + 636/ - 924/ + 552,2 - 96) tt 2 ] /(1296(z - l) 5 ) (3. 



3.2 Counterterms 

The counterterms relevant for calculating on-shell matrix elements are generated by expressing 
the bare parameters in the original Lagrangian in terms of the renormalized quantites. Working 
up to order a s , the only parameters which need renormalization in the present situation are 
the t-quark mass and the b- quark mass (in principle also the s- quark mass if we did not 
work in the limit m s = 0). Using on-shell renormalization for the external 6-quark mass and 
MS renormalization for the (internal) top quark mass, the connection between the bare and 
renormalized masses reads 



m<b,bare 



m t - 5m t , 

777,6 — $ m b , 



5m t 



m t 

5m b a 



a s n 3 
— U F - 
Air e 



47r v e 



6 ln(///m&) + 4 



(3.9) 



Note, these mass shifts not only shift the mass terms like mtit, but also the Yukawa terms 
like ~ g b (rribL — m t R)t where $~ is the unphysical charged Higgs field which appears in 
covariant gauges. These counterterms, induced by the shifts Sm t and Srrib, generate corrections 
for the b — > 57 matrix element, which we denote by 5Mb and SM t , respectively. Writing 
SM f = V5R f (f = t,b) with V given in eq. (POD , we get 



5R b = I ((6,2 - 8) In z - 7/ + 16z - 9) ^ + 4 ln(///m 6 ) + 8/3 
+ (-62 + 8) In 2 z + (20/ - 26z) In z - 19/ + 44z - 25 } 



m w 



2( 



lQ(z-l] 



(3.10) 



while 5Rt is given by 

{6 
- ((18/ + 30/ - 24z) In z - 47/ + 63/ - 9z - 7 ) + 18 z (-3/ -5,2 + 4) In 2 z 

+ (246/ + 114/ - 288z + 96) In ,2 + 44/ - 547/ + 855/ - 413,2 + 61 } x 

(3.11) 



m w 



2, 



24 (z - iy 
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When writing down the expression for SRb, we should mention, that we did not include the 
insertion of Sm^bb in the external 6-quark leg. This is quite in analogy of omitting self-energy 
diagrams for the external legs. Such corrections on the external legs are taken into account 
by multiplying the amputated diagram with the factor \J Z 2t b ~Z%si where Z 2 ^ and Z 2jS are the 
residues taken at the (physical) pole position of the regularized b- and s- quark two point 
functions, respectively. Making use of the expression (in Feynman gauge) 

Z 2 (m) = l-^C F (JL) 

the counterterm 5Mz 2 induced by the Z 2 - factors of the external quark fields reads (again 
writing 6M Z2 = V5Rz 2 ) 

/ a \ 2e ( 2 1 
SR Z2 = - — — (AT 7 oo + e K m ) + - (K 700 + e K im ) 
\m w J Ve ir e 

+ (4 - 6 ln(m 6 /m w ) - ^ lnr) ^ 700 } . (3.13) 



1 2 

- + — + 4 



(3.12) 



4 b — > S7 in the effective theory 

As in the full theory, we first discuss the matrix elements for b — ► 57 of the operators in basis 
(|1.2|) . In section 4.2 we list the various counterterm contributions. 



4.1 Regularized Feynman diagrams 

We write the matrix element M for b — * 57 as a sum of the contributions due to the operators 
0~ in the effective Hamiltonian, i.e., 



i=l 



AiG F X t 

~7T 



Ci(8 7 \Oi\b) 



(4.1) 



To facilitate later the comparison between the results in the two versions of the theory (full vs. 
effective), we write M % = Mq + ^ M\ and cast the term proportional to a s in the form 



a s 
An 



M[ = VRi 



(4.2) 



where V is given in eq. ( 3.2|) . 

We first discuss the contributions of the four- Fermi operators 0\-Oq. As the Wilson coeffi- 
cients of Oi, O3, O4, O5 and 0$ start at order a], we only have to take into account their order 
a® (one-loop) matrix elements; it is well-known that in the NDR scheme only O5 and Oq have a 



non-vanishing one-loop matrix element for b 
0) 



57. Making use of the Wilson coefficients (see 



cm 



47T 



C f 



1 y, 
6 m w 



cm 



47T 



1 , M 3 ~ 
- In — + -E 

2 m w 8 . 



(4.3) 
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R 5 and R 6 are readily obtained 



— — - 



iln^ 

6 m w 



l -E 



Rr 



1 , M 3 ~ 
- In -!— + -E 

2 mvi/ 8 



with 



E 



2, 2 2 15 - 16 z + 4z 2 ), z(18-llz 

— InzH ; - \nz + — -. 

3 6(1 -z) A 12(1 -2 



2 
3 



(4.4) 



(4.5) 



On the other hand, the Wilson coefficient of the operator 2 starts at order a° s . Consequently, 
we have to take in principle one- and two-loop matrix elements of this operator. In prac- 
tice, however, the order a° s (one-loop) matrix element of 2 vanishes and therefore only the 
contribution of the order a\ (two-loop) matrix element remains: 



(4.6) 



Rl 



As this contribution also occurs in the full theory result in section 3.1 (see eqs. ( |3.1|) and ( |3.3| )), 
the explicit expression for the r.h.s. of eq. ( f4.6| ) is not needed for the extraction of C71. 

The order a s contribution of the matrix elements of the dipole operator O7 (see Figs. 4a,b) 
yields 



2<- 



Ri 



3 c 
i° 71 



a 



( —) c 

70 KmwJ In r + In 2 r + 2C 70 In r \n(m b /m w ) -2C 70 \nr . (4.7) 



The first term on the r.h.s. of eq. ( |4.7| ) comes from the tree-level matrix element in Fig 4a, being 
multiplied with the the order a s part (i.e. C71) of the Wilson coefficient C7. The remaining 
terms are due to the one-loop graph in Fig. 4b. Note that C71 is the quantity we ultimately 
wish to extract. Finally, the diagrams of 0$ are depicted in Figs. 4c, d; its contribution is |13| 



Rh 



a 



80 



12 

33 + 2vr 2 + 24 ln(m 6 //i) - 6 i ir 

e 



(4.8) 



4.2 Count erterms 

As the operators mix under renormalization, we have to consider counterterm contributions 
induced by operators of the form Cj SZij Oj. We denote their contributions to b — >• 57 by 

4iG F X t 



SMi. 



V2 



(sjlCJZ^OAb) 



The non-vanishing matrix elements read (using SMij = V SR 



1 .1 1 



5R 



25 



LlfJL 

36 e \m b 



2c 



5R 



26 



2c 



II (JL \ . siirr ,m~ 



4 e \m b 

where we made use of the renormalization constants 



29 1 
27 e 



0:.s 



. $Z 2 5, 5Z 2 s, SZ 2 7, SZ 77 , 5Z 87 ) — — Cp ( — 



1 1 29 4 
12? 4e' 27? e' 



- O70 

e 



4 
3e 



5R. 



87 



(4.9) 

(4.10) 
(4.11) 
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It is well-known that the renormalization of the four-Fermi operators requires the introduction 
of counterterms proportional to evanescent operators f26|| . Calculating b — > 57 up to order 



a s , there are potential counterterm contributions involving evanescent operators needed to 
renormalize 2 . As the initial conditions for the four- Fermi operators (which we partially used 
in section 4.1) depend on the actual choice of the evanescent operators, we have to use the same 
set when calculating their effect of b — > S7. We consistently take both, the initial conditions of 
the four- Fermi operators and the set of evanescent operators from refs. ]T7|, [27], |2"%|. The 



only potentially relevant matrix element of evanescent operators contributing b — > 37 is 

(a7|i^[O a ]|6) , (4.12) 
where the evanescent operator E\ [O2] is of the form 

Ei[0 2 ] = [smlf/SYvlnLc^ c a ^ v ^ u ^Lb ai - (4 + a x e) s ai ^^Lc a2 c a3 Y L K 4 ] K aia2a:jCe4 

■^ol\ol20l^ola 7* ^ a i a 3^ a 2^4 ~^& aia2 ^0304 • (4 .1.3) 

2 

However, as these matrix elements are identically zero (in d dimensions), there are no contri- 
butions from counterterms proportional to evanescent operators. 

Besides the counterterms induced by operator mixing, we also have to renormalize the b- 
quark mass which explicitly appears in the operator 7 and in addition we have to multiply the 



lowest order matrix element by the factor \J Z 2 (m b ) Z 2 (m s ), quite in analogy to the calculation 

in the full theory. The counterterm due to the 6-quark mass renormalization 5M b = V 5Rb 
yields 



SR h 



3 

- + 6 \n(ii/m b ) + 4 
e 



CVo , (4.14) 



when using the on-shell definition for the 6-quark mass, while the counterterm 5Mz 2 = V 5Rz 2 
is given by 

5R Z2 = - [£-) e {^-CVo +- e C 70 + (4-6 1n(m 6 /m w )-^lnr) C 70 } . (4.15) 

5 Extraction of CVi(/iv^t) 

To summarize section 3, the order a s part M[ en of the renormalized matrix element for b — > 37 
in the full theory reads 

M[ en = V \R 1+2 + R 2 + 5R b + SRt + 5R Z2 } , (5.1) 

47T L J 

where the quantities in the bracket on the r.h.s. of eq. ( |5.1| ) are given in eqs. ( |3.5| ), ( |3.3| ), 
(|3.10|) , ( |3.11|) and ( p,13|) , respectively; the prefactor V is given in eq. 



The corresponding renormalized matrix element M[ en in the effective theory can be obtained 
from the information in section 4; M[ en reads 



-=! MT en = V 

47T 



R 2 + R 5 + R e + R 7 + R 8 + SR 25 + SR 2 e + SR 2 j + 5R 77 + 5R 87 + 5R b + 5R Z2 

(5.2) 
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where the various quantities in the bracket are given in eqs. (|4.6|), ( f4.4|) , ( f4.7|) , ( f4.8|) , ( }4.10|) , 
(PH) and fl£Igp . 

Before we are able to correctly extract C71 , a remark concerning the infrared structure is in 
order. We split M[ en into a infrared singular and an infrared finite piece, i.e., 

M[ en = M{f T + M r x %. (5.3) 

As this splitting is not unique (concerning the finite terms), we define the singular part to be 



M^^-(K 700 + eK 701 )^^\nr-2(K 700 + eK 701 )^^ , (5.4) 



where the first and second term on the r.h.s. are due to the two- loop diagrams ( p.5|) and 
the counterterms (|3.13|) , respectively. We do now an analogous splitting for the renormalized 
matrix element in the effective theory, i.e., 

Mr = M[f r +Ml™ n , (5.5) 

with 

M{f r = -C 70 KmwJ lnr-2C 70 KmwJ . (5.6) 

6j r €-i r 

As the matching has to be done in four dimensions, we cannot - strictly speaking - use 
the process b — > 37 to do the matching, because of the infrared singularities. To cancel these 
singularities, we have to include the gluon Bremstrahlung process b — > s^g in both versions 
of the theory. In the effective theory, the process has been worked out in || [Hj, (but the 
explicit result is not important here); the result in the full theory is obtained from the effective 
theory result by replacing C 7 o by K 7Q0 + eK 70 i. The correct physical matching condition consists 
in requiring the infrared finite quantity T = T(b — > 37) + T(b — > s'-fg; E 1 > E™ m ) to be equal in 
both versions of the theory. Due to the specific form of eqs. ( |5.3|) - (|5.6|) and due to the specific 
difference in the bremsstrahlung contribution, it follows that the physical matching condition 
implies 

M{ e J m = MTJ m . (5.7) 
The extraction of C 7 \ is now straightforward. In summary: Writing the Wilson coefficient 
C 7 ([iwt) at the matching scale fiwt in the form 

C 7 (fiwt) = C 70 ([iwt) + ~r~ C 7 i([iwt) , (5.8) 

we obtain (in the naive dimensional regularization scheme) 

, 2z(8z 3 + 61z 2 -40z + 4) T . 1 2z 2 (3z 2 + 23z - 14) , 2 

= 9(*-l)< - Ll(1 "I )+ 3(z-l)« 

2 (51z 5 + 294z 4 + 11582 3 - 16972 2 + 7422 - 116) 



81(z - l) 5 

15202 4 + 129612 3 - 121262 2 + 3409^ - 580 



In z 



486(2- l) 4 
Az 2 (3z 2 + 23z - 14) 



3(2- l) 5 

2 (106z 4 + 2872 3 + 1230z 2 - 12072 + 232) 



In z \n(jj wt /mw) 

\n(/j, wt /m w ) . (5.9) 
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Here, z = (m t (^wt) I ' m w) 2 > where m t (fj, wt ) is the MS top quark mass at the renormalization 
scale nwt- The lowest order function CV is given in eq. (|2.24|) . 

Taking into account that the result of Adel and Yao [ffl is given in the so-called R* renor- 
malization scheme, we got the same result for Cn(^wt)- 



6 b — > sg in the full theory 

As in the h — > S7 case we first give the results for the two-loop diagrams and then move to the 
counterterm contributions. 



6.1 Two- loop Feynman diagrams 

We again write the b —>■ sg matrix element M in the form M = Mq + ^ Mi. Using the 
"reduction technique" described in section 2.3, the complete set of two-loop Feynman graphs 
is given by the abelian diagrams in Fig. 2 and by the non-abelian diagrams in Fig. 3, which 
involve the triple gluon coupling. The crosses in Fig. 2 and Fig. 3 show the possible locations 
from where the gluon can be emitted. Of course the graphs with a cross at the W line in Fig. 
2 have to be omitted. Working in the limit m s = 0, diagram 5b in Fig. 2 vanishes. It is 
convenient to write Mi in the form 



O'.s 



Mi = W Ql+ 2 - Ql+ 2 - Ql , (6.1) 



47T 

where the quantity W is defined as 

W = —j=— — ( S1 \0 8 \b} tree . (6.2) 

In eq. ( |6.1| ) <5} +2 denotes the sum of the first and second contribution in the Heavy Mass 
Expansion for an internal quark of flavor / (/ = t, c); Q z c is the third contribution in this 
expansion, which only has to be considered for the light internal quarks. Like E? c in eq. ( |3.1| ) 
of section 3.1, Q\ is just 

Qc = "4 , (6.3) 
where Q2 is the quantity defined through the relation 

(sg\0 2 \b) = ^(sg\0 8 \b) tree Q 2 . (6.4) 

As exactly the same term also appears in the effective theory, Ql drops out when extracting 
the 0(a s ) correction to the Wilson coefficent C 8 . 

The dimensionally regularized expressions for Q 1+2 = Q\ +2 — Ql +2 can be written in the 
form 

Q 1+2 = UK 800 + eK 80 i)^^\nr-3(K 800 + eK 80 i + e 2 Kso2)^- 

/ \ 2e 

3 ( — ) 

-- (K 800 + eK 80 i) [2 + lnr - A\n(m b /m w ) + 2m] 
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, lc 



+hi w h h 2 In 2 r + h 3 In r ln(mb/miy) + h± lnr 



e 

,2 



+/i 5 ln(mb/mvi/) + /i6 In (m 6 /mvy) + . (6.5) 



The first term on the r.h.s. of eq. (6.5) is due to infrared singularities coming from the 
(abelian) graph in set 3 in Fig. 2, where the gluon is radiated from the internal quark; the 
infrared structures appearing in the second and third term are due to non-abelian diagrams in 
Fig. 3. Eq. ( J6.5| ) shows that the infrared singularities again just multiply the dimensionally 
regularized version of the lowest order matrix element (see eqs. ( |2.20| ) and ( [2.21|) ). 



The functions K 800 and K 80 i appearing in eq. ( |6.5|) are given in eqs. (|2.25|) and ( p.27| ) . We 



note that the function K 802 is not needed explicitly in order to extract C 81 , as we will see later. 
The functions hi in eq. (|6.5| ) read (z = {m t / mw) 2 ) 

z (774/ + 810z + 144) In z + 137/ - 823/ + 257/ - 425/ + 958,2 - 104 23 
1 72 (z - l) 5 27 

(6.6) 

2 8 8 

^2 = 2 ^soo , h = -C K 800 , h 4 = -- K m0 , h 6 = -6 K 800 (6.7) 

z (162z - 72) In z + 11/ - 110/ + 57/ + 82^-40 „ . T ^ 
K = _ 1 + 6t7cK 80Q -2 (6.8) 

z(Az 3 -40/ -Alz- l)Li(l - ±) 8. T ^ 59 2 185 

n>7 = : rz ~ % 7T K 8m 7T 

7 6(z-l) 4 3 108 324 

- (359642 s + 54756/ + 2592z) In 2 z 

+ (7452/ - 42660/ - 92772/ - 73164/ + 48984z - 3360 + 3186 tt 2 (/ - /)) In z 
+ (844/ + 40012/ + 90580/ - 148588/ + 16688z + 464) 



+ (-885/ + 3363/ - 9381/ + 13629/ - 8142^ + 1416) ir 2 /(2592(z - 1 



(6.9) 



6.2 Counterterms 



As the discussion concerning the counterterms induced by the shifts in the t- and b- quark 
masses is exactly the same as in the b — > S7 process in section 3.2, we give immediately the 



result. Writing 8Mf = W5Qb (f = t, b) with W given in eq. (|6.2|) , we get 



SQ b = {- (2 lnz + z 2 -4z + 3) + 4 ln(/x/m b ) + 8/3 ) 

+2 In 2 z + 2z (z -4) lnz - z 2 + 8z - 7} f— V* n , - 1X , (6.10) 
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6 

5Q t = \- (-6z(z + l) Inz + z 3 + 9z 2 -9z-l) + 18z{z + l) In 2 z 



+ (-6z 3 - 84z 2 - 18z + 24) In z + 5z 4 - 10z 3 + 126z 2 - 158^ + 37 } x 

(6.11) 



\m w J 3 (z — l) 5 

Also the counterterms due to the \fZ% factors of the external quark fields are obtained in the 
same manner as in section 3.2, leading to (5M Z2 = W 8Qz 2 ) 

oQz 2 



mw 



) - { — (K 800 + e K m ) + - (K 800 + e K ml ) 
/ 3 Uir e 



+ (4 - 6 \n{m b /m w ) - - In r \ K 800 | . (6.12) 

For the b —>■ sg case there are additional counterterm contributions due to the strong 
coupling constant renomalization and due to the \fZ~z factor associated with the external gluon. 
Denoting the combined effect by 5M g = W 5Q g , one obtains 

8Q 9 = (~- + f) (^soo + e^soi) • (6.13) 



As the finite term / will appear also in the corresponding counterterm in the effective theory, 
it will drop out when extracting Csi- 



7 b —> sg in the effective theory 

7.1 Regularized Feynman diagrams 

In the effective theory the matrix element M for b — > sg is of the form 

8 4?T y A 

M = £ikT ; M* = -^Ic^sgp^b) . (7.1) 
i=\ V2 

We write M % = Mq + ^ M\ and put the term proportional to a s into the form 

-JLM\ = WQ % , (7.2) 

where W is given in eq. (|6.2j). As the discussion how to get the quantities Qi is basically 
identical as in the b — > sj case in section 4.1, we just give the results. Among the four- Fermi 
operators, only O2 and O5 yield non- vanishing matrix elements for b —* sg. We get 

Q2 , Q 5 = -\\n-^-- l -E , (7.3) 
9 m w 6 

where E is given in eq. ( |4.5|) . Again, we do not have to know Q 2 explicitly, because this term 
also appears in the full theory result; it drops out when extracting Cgi- 
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While there is no contribution from the dipole operator O7, there are various diagrams 
associated with the operator 8 (see Figs. 5,6). The sum of all these contribution is given by 



1 3 
7; ^80 mr — 6 Oso 2 9 ° 80 



+Q 



so 



(^) 2£ 8 2 

— V m Hv 1_ g 2 77 In (m&/mw) i tt H — In 2 r — 6 \v?(mb/mw) 

3 6 3 3 



-f lnr + I lnr ln(m ft /m w ) - ~ ln(m b /m w ) + ~ + ^ vr 2 ) + C 81 . (7.4) 

3 3 6 6 6b J 

When comparing with the full-theory expression Q 1+2 in eq. (|6.5|) , one immediately realizes 
the correspondence of the infrared singularities. To this end it is important that one carefully 
disentangles everywhere infrared and ultraviolet poles. Especially, one should use the formula 

d d r 1 i n 1 \ . . . f d d r 1 , . 

I instead of / - — -7 = . (7.5) 



(2n) d {r 2 ) 2 16tt 2 Ve e ir J J {2Ti) d (r 2 ^ 2 

An example, where such a situation occurs, is the diagram in Fig 6c. 



7.2 Counterterms 

As the operators mix under renormalization we have to consider counterterm contributions 
induced by operators of the form 5Zy Oj. We denote their contributions to b — > sg by 

5M tJ = A%( ^ Xt (sgldSZijOjlb) . (7.6) 
The non- vanishing matrix elements read (using SMij = W 5Qij) 

0Q25 = - n - > ^V28 = ^ - > d V88 = ^~ ~ <~-80 , (7-7) 

9 e \m b ) 27 e 3 e 

where we made use of the renormalization constants [|J 

9e 47T 27e An 3e An 

We note that there are no contributions to M(b — > sg) from counterterms proportional to 
evanescent operators. 

In analogy to the 6 — > 57 case in section 4.2, there are the counterterms from renormalizing 
the 6-quark mass which explicitly appear in the definition of the operator 8 and from the 
\[Z~2 factors for the external quarks. The counterterm due to the 6-quark mass renormalization 
5M b = W8Q b yields 

4 T3 

$Qb = -r - + 6 ]n(n/m b ) +4 
o L 6 



C 80 , (7.9) 
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when using the on-shell definition for the 5-quark mass ( |3.9|) , while the counterterm 5Mz 2 = 
W SQz 2 is given by 

5Qz 2 = - (—T \ {— C 80 + - C w + (a - 6 hx{m b /m w ) - \ lnr ) C 80 ) . (7.10) 

Finally , there are counterterms due to the strong coupling constant renormalization and due 
to the \fZl of the external gluon. As in the full theory, we only give the combined counterterm 
6M g = WQ 9 

5Q 9 =(-- e +f)c 80 ■ (7.11) 

As / is the same finite quantity as in the corresponding result (|6.13 ) obtained in the full theory, 
we do not need its explicit form, because it drops out when extracting C%\. 



8 Extraction of Csid^wt) 



To summarize section 6, the order a s part M[ en of the renormalized matrix element for b — > sg 
in the full theory is given by 

?L M[ en = W \Q 1+2 + Q 2 + 5Q b + 6Q t + 5Q Z2 +6Q g ] , (8.1) 
Att 1 J 



where the quantities in the bracket on the r.h.s. of eq. (|8.1| ) are given in eqs. ( |6.5| ), (|6.3|) , 
(|67TOp , (|6U1|) , (|67ll ) and (gig) respectively; the prefactor W is given m eq. ( p.Z\ ). 

The corresponding renormalized matrix element in the effective theory can be obtained from 
the information in section 7; M\ en reads 

rv ~ r s\ ~ ~ ^ ^ ~ 1 

-H- M[ en = W Q 2 + Q 5 + Q 8 + SQ 25 + 5Q 28 + 5Q 88 + 5Q b + 5Q Z2 + 5Q g , (8.2) 
where the various quantities in the bracket are given in eqs. ( |7.3|) , ( |7.4| ), ( |7.7| ), (|7.9|) , ( |7.10| ) and 

Before we extract C 8 i, which enters M[ en via Q 8 (see eq. ( |7.4|) ), we should point out that 
the discussion concerning the infrared singularities is similar as in the b — > 57 case in section 5; 
all the formulae are written in such a way that we simply can discard the terms proportional 
to the poles in ej r in both versions of the theory. The extraction of C 8 i is then straightforward. 

To summarize: Writing the Wilson coefficient C 8 (fiwt) at the matching scale \iwt in the 
form 

01 

C%{^wt) = C 80 (fiwt) + -7— Csi(nwt) , (8-3) 

we obtain (in the naive dimensional regularization scheme) 

, z(Az 3 -AOz 2 -Alz-1) T . n 1 z 2 (17z + 31) l2 

° M = Ll(1 -^- 2(*-l)» 

210z 5 - IO862 4 - 4839z 3 - 3007z 2 + 2UAz - 304 



216(2- l) 5 
6II2 4 - 13346^ 3 - 295952 2 + 1510z - 652 
+ 1296(2 - l) 4 + 
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In z 



z 2 (17 'z + 31) , , u wt 89z 4 - U6z 3 - U37z 2 - 95(k + 152 , u wt 
H ; 77 — In 2 m 1 ; — m . 

(-2-1) 5 m w 54(z-l) 4 m w 

(8.4) 

Here, z = (mt(fxwt) / '™>w) 2 , where m t (fiwt) is the MS top quark mass at the renormalization 
scale fiwt- The lowest order function C$o is given in eq. ( |2.25| ). 

Taking into account that the result of Adel and Yao [Q] is given in the so-called R* renor- 
malization scheme, our result is identical. 



9 Summary 

The order a s corrections to the Wilson coefficients CV and C% are a very crucial ingredient for 
the prediction of the branching ratio for b —>■ Xgj in next-to-leading logarithmic precision. As 
these corrections, which involve many two-loop diagrams in the full theory, have been calculated 
so far by one group [l[] only, we presented in this work a detailed recalculation. We extracted 
the 0(a s ) corrections to C-j and C§ by comparing the on-shell processes b — > 57 and b —>■ sg in 
both versions of the theory. We evaluated the two-loop integrals in the full theory by using the 
Heavy Mass expansion method. Our a s corrections (C71 and Csi) to the Wilson coefficients C-j 
and Cs completely agree with the findings of Adel and Yao. 

We should point out that our result (as well as Adel and Yao's) for C n (fj,wt) and C 8 i(fiwt) 
is a priori specific to the basis given in eq. ( |1.2j) . However, the same answer is obtained for 
these Wilson coefficients when working in the basis recently used by Chetyrkin, Misiak and 
Miinz. 
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Figure Captions 
Figure 1 

a Lowest order diagram for b — > 57 in the full-theory. A cross denotes a possible location where 
the photon can be emitted. The wavy line stands for a W or unphysical Higgs boson ($). 
In the b — > sg case the cross at the W 7 / $ has to be ignored. 

b Typical two-loop graph for b — > 57. 

c-e Subdiagrams of b) which contribute in the Heavy Mass Expansion. See text. 



Figure 2 

Complete list of two-loop diagrams for b — > 37 in the full-theory. A cross corresponds to a 
possible location for the photon emission. 

For the b — > sg process, this figure is a complete list of diagrams not involving the gluon triple 
coupling. (In the b — > s# case the crosses at the wavy (W/<&) lines should be ignored.) 

Figure 3 

Complete list of two- loop diagrams involving the triple gluon vertex (for the b — > sg process). 
Figure 4 

Diagrams associated with the operators 7 and 8 in the effective theory for b — > 57. See text. 
Figure 5 

Abelian diagrams associated with the operator 0$ in the effective theory for b — > sg. See text. 
Figure 6 

Non- Abelian diagrams associated with the operator 0$ in the effective theory for b ^ sg. See 
text. 



24 





Figure 1 




Figure 2 




Figure 3 




Figure 4 




Figure 5 




Figure 6 



